Theory of cooling by flow through narrow pores 
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We consider the possibility of adding a stage to a dilution refrigerator to provide additional 
cooling by "filtering out" hot atoms. Three methods are considered: 1) Effusion, where holes having 
diameters larger than a mean-free path allow atoms to pass through easily; 2) Particle waveguide- 
like motion using very narrow channels that greatly restrict the quantum states of the atoms in 
a channel. 3) Wall-limited diffusion through channels, in which the wall scattering is disordered 
so that local density equilibrium is established in a channel. We assume that channel dimension 
are smaller than the mean-free path for atom-atom interactions. The particle waveguide and the 
wall-limited diffusion methods using channels on order of the de Broglie wavelength give cooling. 
Recent advances in nano-filters give this method some hope of being practical. 



We investigate here the possibility cooling a gas by "fil- 
tering out" hot atoms, perhaps as an adjunct to a dilution 
refrigerator with solutions of 3 He in liquid He.[lJ Our 
method involves passing a degenerate Fermi gas through 
narrow constrictions formed by pores in a partition. Un- 
der certain conditions we find that this particle "leakage" 
allows a lowering of the temperature of the remaining gas. 
Our study adds to many recent experimental and theo- 
retical analyses of quantum size effects in the behavior of 
particles in nanopores (as in Refs. 

An initial idea of how to remove the hot atoms is 
suggested by a technique used in electron hetero junc- 
tion physics, 0, 0, Q, where the electron gas is passed 
through a narrow constriction formed by a gate potential. 
Because the constriction is narrow the bands of states al- 
lowed in this "particle waveguide" are widely separated, 
which means that not all energies are allowed through. 
Adjusting the states in the channel can allow selective 
passage particles in states at the Fermi energy, so that 
one removes only hot atoms. In such an approach the 
constriction must be of order of the de Broglie wave- 
length of atoms at the top of the Fermi surface, which is 
roughly the separation between fermions. Because of the 
limited states (that is, bands) we characterize channels 
of this size as "narrow." We will see that this approach 
can be made to work under appropriate conditions on the 
nano-pores. 

An alternative possibility involves channels or pores 
with diameters much larger than the de Broglie wave- 
length but still less than the mean-free path of the 
fermions in the gas, which can be large due to its 1/T 2 
behavior. We call these "wide" pores. 

We will consider different kinds of flow through the 
holes of the two sizes mentioned: effusion, waveguide 
flow, and wall-limited diffusion. In effusion, the holes 
are by definition wide and the states in the channel re- 
main three-dimensional with no banding. Moreover the 
walls are sufficienly smooth that the particles undergo no 



back scattering nor do they come into equilibrium with 
the channel walls. Effectively all that the pores do is to 
allow particles already directed in the positive z-direction 
to pass through the membrane. We might hope that this 
would allow cooling because the intensity of fast atoms 
passing through a hole is larger than that of the slow 
atoms and these, on average, carry more energy. Indeed 
it is a standard textbook exercise^^ to show that this 
works for a Boltzmann gas. But does it work for a de- 
generate Fermi gas? 

In waveguide flow the channels are narrow enough to 
have well-defined bands. Again the particles are assumed 
to undergo no back scattering nor do they come into ther- 
mal equilibrium with the channel walls. Nevertheless be- 
cause of the bands only particles with certain energies 
are allowed in the channel. If, for example, the Fermi 
energy of the gas in the container is coincident with the 
bottom of the lowest band one might expect that only 
high energy particles would get through the holes and 
the remaining gas would be cooled. 

A third situation is wall-limited diffusion or Knudsen 
flow, which could occur in wide or narrow pores. There 
is scattering at the walls including back scattering. The 
rate of diffusion in the channel depends on the diam- 
eter of the pore, the particle velocity, and the density 
gradient that maintains the flow. In this and the above 
methods we will assume the density difference is main- 
tained by pumping away the particles that pass through 
the membrane. Both effusion and wall-limited diffusion 
were considered by one of the authors some time ago for 
enhancing polarization. 

If the Fermi gas in the cooling cell is, say, a 1% solution 
of 3 He in liquid 4 He at millikelvin temperature, then the 
Fermi temperature is 124 mK, the de Broglie wavelength 
about 3 nm, and the mean-free path is about 1 /jm at 
T = 15 mK. 

Suppose the number of atoms in a box (Bl) is N with 
density n. The atoms pass through a membrane and en- 
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ter a second box (B2) maintained (by pumping) at a 
much smaller particle number N% and density n^. The 
membrane contains a great number M of holes each of 
diameter d and of total hole area A. (See Fig. 1.) We 
will compute particle current Jn and energy current Je 
passing through the holes. The rates of particle and en- 
ergy change, dN/dt = AJm and dE/dt = AJe in the 
container Bl combine in 
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to give the cooling power CydT/dt where Cy is the heat 
capacity at constant volume in Bl. 
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FIG. 1: Schematic diagram of the apparatus. 3 He in solu- 
tion with liquid 4 He enters the cooling chamber Bl, passes 
selectively through the membrane into B2. A gas at lowered 
temperature remains in Bl. The 3 He density in B2 is kept 
low by pumping. The gas is recycled, after being recooled by 
a dilution fridge, back into Bl. 

In the case of effusion the particle current is 
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where n p is the momentum distribution function and 
the factor of 2 accounts for spin degeneracy. There is 
a similar expression for Je- For a classical gas passing 
through wide holes we find a cooling power given by 

1/2 

— nfcsAT {HbT /2'Km) . The minus sign implies cool- 
ing in agreement with the usual textbook treatment. [To| 
In the degenerate limit effusive cooling fails. 

For particle waveguide motion (in narrow holes) we 
need to take into account the banding of the states due to 
the limited transverse motion; the energies are given by 
e n2 = e n +e z with e z — h 2 k 2 /2m. We will assume a square 
cross section with width d so that e„ = eor/ n where eo = 
TT 2 h 2 /(md 2 ) and rj n = 1, 2.5, ...(n 2 , + n 2 )/2,.... (Using a 
circular cross section gives no significant differences in 
our final results.) In this case Eq. J5J is replaced by 
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with n the chemical potential (w ee at low T). In 
terms of the one-dimensional Fermi integrals Fi(7) = 



r o °° dzz l l/(e z -~* + 1), the cooling power for the waveg- 
uide case is 



where 
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The filling of the bands in the channels is determined by 
a n = /3(/t-e n ) with /3 = (fc s T)- 1 and £ = {dE/dN) T V 
is a property of the gas in the cooling chamber (Bl). For 
cooling pw must be positve. If 7 = P/j, and t = T/Tf 
then the relation between density and chemical potential 
leads to the usual transcendental equation for fi : 



^3/2 Fi/2(7) = j 

Further, taking the derivative of E versus N gives 



£3/2(7) 
£1/2(7) 
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where G t (i) = J °° dzz l e z ~^ /(e z ~~< + l) 2 = IFi-i. Given 
a value of t from the temperature and the concentration, 
we must solve Eq. © for 7; we put that into Eq. JJJ to 
get £. We find it convenient to introduce the the ratio 
of lowest band edge to Fermi energy y = £q/£f so that 
P^ = yr/n/t. 

For numerical calculations we fix T at 15 mK as a rea- 
sonable incoming temperature in the cooling cycle. Then, 
if the number of holes in the membrane is M =10 , the 
prefactor in Eq. (0J is 13 ynW. In the case of wide holes 
(where waveguide flow becomes effusion), the sum over 
bands in Eq. (@J becomes a double integral over trans- 
verse momenta and in this continuum limit we find 

7T t 

pw ~ [^2(7) — £-^1(7)] continuum limit (8) 

2 y 

We have made no assumptions here about whether the 
system is degenerate or not — just that the channel states 
are now continuous. If the factor in square brackets is 
positive, then it says that the smaller y, the larger is the 
cooling. This result stems from just having bigger and 
bigger holes, allowing more hot atoms out. But of course 
there is a limit to how big the holes can be to maintain 
a pressure differential across the membrane. 

Next consider the highly degenerate limit of the last 
form, for which F n — 7™/n and £ = 7, so that 
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7 3 degenerate continuum limit (9) 



which agrees with a direct calculation — no cooling occurs 
in that case, because there are many channel states below 
the Fermi energy allowing low energy atoms to escape. 
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Indeed the gas left behind can then end with a higher 
temperature. 

For the fully classical limit we take e 7 to be very large 
so that 
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classical continuum limit. 
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(10) 



This result again agrees with the textbook classical effu- 
sion result. 

In Fig. 2 we show numerical results based on Eq. It 
is useful to plot ypw rather than pyy because then all the 
plots collapse onto a single curve for t > 1, as implied by 
Eq. iJSJ. The t dependence for large t is t~ x / 2 and we are 
indeed in the asymptotic region described by Eq. (|10fl . 
We see that y = 1 is a special point with the lowest band 
edge lined up with the Fermi energy. (Recall however, 
that the cooling power prefactor, which is not included 
in pw contains a T 2 so all cooling powers go to zero with 
decreasing T.) When the lowest band edge is lower than 
the Fermi energy (y < 1) we lose cooling at very low T 
(as in Eq. (0), but surprisingly we regain it for larger 
T due to the contributions of the higher bands . Finally 
for all band edges above the Fermi energy (y > 1) we get 
cooling for all T, but, because the curves in the figure 
are multiplied by the factor y, the actual cooling power 
is diminished when we divide out this factor, and these 
cases are not as useful. 
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FIG. 2: Reduced waveguide cooling power as a function of 
temperature parameter t for various band-edge y values for 
a square cross section. Because of the multiplicative factor y 
all curves collapse onto the same classical limit curve at large 
t. 



The parameters we have used are not independent. Be- 
cause Tp depends on a; 2 / 3 , where x is the concentration, 
and eo depends on d~ 2 , then for T = 15 mK we find the 
relation 



t = 1.8 x 10~ 2 y (d(nm)) 2 



(11) 



From the plot, we see that if we have y < 1 we get a 1/y 
enhancement in cooling power, but to gain that we must 
also have t > 0.4 to avoid heating. For y > 1 we can 



go to lower t values but we lose cooling power because of 
the 1/y factor and because of the dip in the curves. 

Consider d = 10 nm, M — 10 11 which is perhaps 
within practical reach. The value t — 0.4 (Tp — 35 mK 
and x — 0.2%) implies y — 0.2 and the cooling would be 
on the order of a 50 /iW in this ballistic-flow waveguide 
case. Such a result would be quite a remarkable cool- 
ing rate, but the assumption of perfectly smooth walls is 
optimistic. However, the inevitable coating of the pore 
walls by a couple of layers of solid 4 He will enhance the 
smoothness. 

It perhaps seems more likely that the walls of the chan- 
nel would cause scattering, including back scattering, in- 
tcrband transitions, etc. Thus we consider next a sim- 
ple wall-limited diffusion or Knudsen flow model. Our 
starting point is a kinetic equation for the distribution 
function in collision-time approximation^!] >E2l 
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where r is the time between collisions with the wall. The 
left side can be written as — v pz dn p ^ /de p dfj,/dz, where 
is the local equilibrium distribution function and 
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8n p is the correction to local equilibrium. The gradi- 
ent in chemical potential dp/dz is proportional to the 
gradient in density. A much more rigorous approach to 
such a kinetic equation is described in Ref. [l3j . We con- 
sider again the case of a very narrow channel containing 
banding. We solve for 5n p and use it to compute, say, 
the particle flux as 
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From Eq. Q the cooling power is found to be 
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where we have introduced the same prefactor as in the 
waveguide case so that 
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with Cn = — yVn/t- In this equation the collision time r 
is writen in terms of the transverse velocities in a channel. 
That is, we write r ~ d/(2v) with v = y/2 (e n ) jm where 
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(16) 

If we divide this quantity by the lowest band edge, then 
we have (e rl ) /eo = (Vn) used in Eq. (|15|l . 

We have to evaluate the derivative d-f/dz. What we 
mean by this quantity is dj/dz = (dj / dn) (dn / dz) since 
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it is the gradient in density n that drives the flow. We 
can find this by taking the derivative of the self-consistent 
expression, Eq. JBJ. In that equation t depends on n 
because ep does. We find 
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where in the middle form we have taken the density gra- 
dient as n/L with L the length of a channel. What we 
compute then is 
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(18) 

We again plot the product, y pk- The results are 
shown in Fig. 3; note that the ratio d/L is not included 
in the plots. The curves differ from those of the waveg- 
uide in a several ways. The curves are always positive, 
i.e., represent cooling, and the y = 1 degenerate case 
diverges (the T 2 in the prefactor cures that). We have 
checked the numerical calculations by doing various lim- 
iting relations analytically. For example, the continuum 
(small y) case in the limit t — * is pk — > 2Ald/Ly. 



In summary, we have examined here the physics of par- 
ticle flow through narrow pores and estimated the pos- 
sibility of cooling by this method and have found some 
potential for success. As nano-technology improves, the 
possibilities may increase. A very smooth-walled chan- 
nel that would provide waveguide type flow would give 
the greatest cooling. The more probable situation of 
Knudscn flow, while providing cooling over all parameter 
ranges, has the factor d/L reducing the cooling power. 
However, even that circumstance does not make it im- 
possible. The numerical results give hope that this ap- 
proach can lead to a add-on device to extend the range 
of a dilution refrigerator. Experiments are being planned 
to test the potential of the method. 14] While we have 
considered the possibility of a practical application of 
nanopores here, experiments on this kind also provide 
interesting physics, namely detecting quantum size ef- 
fects in the narrow channels and the resultant restriction 
of states as already seen experimentally. 0, El 
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FIG. 3: Reduced Knudsen cooling power as a function of 
temperature parameter t for various y values for a square 
cross section. The factor of d/L in the overall cooling power 
is not included here. 

As in the waveguide case, small y gives amplified 
cooling. Indeed, for M — I0 11 holes, y = 0.1 gives 
PxL/d « I.6/2/ — 16. Again consider d = 10 nm. The 
cooling power is 13 /iWxl6d/L = 208g?/L ^iW. A mem- 
brane width of L = 1 /im gives a cooling of just 2.1 /LW. 
By Eq. I|ll(l we have t = 0.18 with the y value chosen or 
Tp = 83 mK and x — 0.5%. This cooling power might 
be useful if the width L used is not too optimistic. The 
cooling would be enhanced if a larger M value were avail- 
able. For small y values the continuous degenerate limit 
gives an upper limit on the cooling power within these 
conditions. We have 
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